|
In numerical linear algebra, the Jacobi eigenvalue algorithm is an iterative method for the calculation of the eigenvalues and eigenvectors of a real symmetric matrix (a process known as diagonalization). It is named after Carl Gustav Jacob Jacobi, who first proposed the method in 1846, but only became widely used in the 1950s with the advent of computers. == Description == Let ''S'' be a symmetric matrix, and ''G'' = ''G''(''i'',''j'',''θ'') be a Givens rotation matrix. Then: : is symmetric and similar to ''S''. Furthermore, ''S′'' has entries: : where ''s'' = sin(''θ'') and ''c'' = cos(''θ''). Since ''G'' is orthogonal, ''S'' and ''S''′ have the same Frobenius norm ||·||F (the square-root sum of squares of all components), however we can choose ''θ'' such that ''S''′''ij'' = 0, in which case ''S''′ has a larger sum of squares on the diagonal: : Set this equal to 0, and rearrange: : if : In order to optimize this effect, ''S''''ij'' should be the off-diagonal component with the largest absolute value, called the ''pivot''. The Jacobi eigenvalue method repeatedly performs rotations until the matrix becomes almost diagonal. Then the elements in the diagonal are approximations of the (real) eigenvalues of ''S''. 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「Jacobi eigenvalue algorithm」の詳細全文を読む スポンサード リンク
|